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Abstract 

Recent advancement in the molecular Rayleigh scattering based technique allowed for simultaneous 
measurement of velocity and density fluctuations with high sampling rates. The technique was used to 
investigate unheated high subsonic and supersonic fully expanded free jets in the Mach number range of 
0.8 to 1.8. The difference between the Favre averaged and Reynolds averaged axial velocity and axial 
component of the turbulent kinetic energy is found to be small. Estimates based on the Morkovin’s 
‘Strong Reynolds Analogy’ were found to provide lower values of turbulent density fluctuations than the 
measured data. 

I. Introduction 

Reynolds averaging involves dividing a fluid flow parameter w into a time-averaged w , and a 
fluctuating w' part: 

1 T 

w (t) = w + w'(t), w=Lim— wdt (1) 

T— >oo T J 
0 

In the Favre-averaging process a density-weighted filtering is used to calculate time average w and the 
remaining fluctuating part w": 


w (t) = w + w "(t) , 


w = — Lim — 

p T— >oo T 


T 

Jp w dt 
0 


( 2 ) 


where p denotes air density. As there is no mass flux across the Favre-averaged streamline, the use of the 
Favre-averaged variables makes the governing equations for the mean density, mean velocity and the 
mean enthalpy more compact. In compressible flow, such simplifications are unattainable by the 
Reynolds averaging process. Experimental data, such as those measured by Particle image Velocimetry 
and the present Rayleigh scattering technique, however, produce Reynolds averaged quantities. The 
differences between the two averages for parameters such as the axial velocity u and the turbulent kinetic 
energy become a concern whenever computational results are to be compared against experimental data. 


NASA/TM— 2005-213564 


1 



The present effort is to determine the differences experimentally by measuring correlations between 
density and axial velocity fluctuations. 

The measurements have been possible due to the recent advances in a molecular Rayleigh scattering 
based technique to simultaneously measure density and velocity fluctuations in high-speed flows 
(Seasholtz, Panda, and Elam, 2002). Fluctuations occurring over a frequency range of 0 to 50 kHz have 
been measured. The technique is based on laser light scattering from the gas molecules present in air. Since 
neither any seed particles nor any intrusive probes are used, the technique is strictly non-invasive. 
Implementation of the technique requires special attention to the cleanliness of the air stream, isolation of the 
sensitive optical components from jet noise, usage of a Fabry-Perot interferometer and low-level light 
measuring electronics. Description of the technique is provided later in the text. 

The present work is part of a database development project and the free jets used for the current 
investigation were also the subject of some of the earlier studies. The time averaged velocity u , 
temperature T , and density p surveys from these jets have been reported by Panda, Seasholtz, Elam, and 
Mielke (2004). Schlieren photographs and various statistics of density fluctuations were provided in 
Panda and Seasholtz (2002). The far field noise sources were identified by correlating turbulent density 
and velocity fluctuations with far field noise pressure fluctuations in Panda, Seasholtz, and Elam (2003) 
and Panda and Seasholtz (2002). 

A second goal of the present work is to validate Morkovin’s strong Reynolds analogy in unheated 
jets. The origin of Favre averaging to remove density fluctuations from the time-averaged equations of 
motion lies with this hypothesis. A physically intuitive description of the hypothesis can be found in a 
review article by Peter Bradshaw {1911): “ It is generally accepted that the direct effects of density 
fluctuations on turbulence are small if the root-mean-square density fluctuation is small compared with 
the absolute density: this is Morkovin ’s hypothesis . . . This means that the turbulence structure of 
boundary > layers and wakes at free-stream Mach numbers . . . less than about 5, and of jets at Mach 
numbers less than about 1.5, is closely the same as in the corresponding constant-density > flow.” 

“In compressible turbulent flows the velocity >, temperature and density... all fluctuate; in analytical 
work it is more convenient . . . to treat the equivalent fluctuations of vorticity, entropy (or total 
temperature), and acoustic pressure. Morkovin ’s important contribution . . . was to use the limited data 
available in 1961 to show that in non-hypers onic boundary layers the acoustic mode is negligible. ..and 
the entropy (total temperature) mode is very > small for conventional rates of heat transfer ... It follows 
from Morkovin ’s findings that 


?r«l, 

P 


fo_ 

To 


«1 


T' 


so that (y-l)M 2 

p T u 


( 3 ) 


The direct quote ends at the equation. In the above, p is static pressure, T 0 is total temperature, p is gas 
density, T is static temperature, y is the ratio of specific heats, u is the axial component of velocity and M 
is the local Mach number based on the local velocity and the local speed of sound. 

Validation of the Morkovin’s hypothesis and measurements of the Favre averaged variables has been 
attempted in the past using hot-wire anemometry. The signal analysis from hot-wires becomes involved 
and carries various assumptions for applications in compressible flows. The advancement of the optical 
measurement techniques brought new tools better suited for the high velocity flows. Bowersox (1996) 
used Faser Doppler Anemometry in conjunction with hot-wires to provide a more straight forward 
measurement of Favre averaged velocity and Reynolds stresses. These efforts were mostly confined to the 
wall-bounded shear layers. In free shear layer hot-wires break easily and the authors do not know of any 
earlier study addressing the present goals. 
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Difference Between Favre- Averaged and Reynolds-Averaged Axial 
Velocity and Axial Component of Kinetic Energy 

Following the definitions of Reynolds averaging and Favre averaging shown in equations (1) and (2) 

g-PU- fr'+pXu' + iQ -u , (4) 

P p P 

Flence, u-u = ^- (5) 

P 

Therefore, the differences in the time averaged velocity is due to the correlations between the density and 
velocity fluctuations. A similar equation holds for the fluctuating terms. Since u+u"=u+u', 


u — u = 


pu 


( 6 ) 


Next kinetic energy is considered which involves product of fluctuations. For example Reynolds- 
averaged kinetic energy k = 0.5(u'u'+ v'v'+ w'w'), where v and w represent respectively, radial and 
tangential components of velocity. In the present program only the axial component of velocity was 
measured. The axial component of kinetic energy in Favre- averaged fluctuations: 


Jy _ pu V _ (p + p')uV _ | p'uV 


( 7 ) 


Using equation (6) the first term on the right hand side can be expressed as: 


uV=u'u'+ 


^p'u '^ 1 


V 1" J 


The second term on the right hand side of equation (7) can also be written as: 


p u u _ p uu 


-2 


pu 


V r J 


( 8 ) 


( 9 ) 


Flence, the difference in axial component of kinetic energy: 


2 


uV'-u'u 

V J 


p'u'u' 


f P'u'V 


( 10 ) 


Above equation shows that in addition to density-velocity correlations, density-velocity 2 correlations have 
significant contributions to the difference. In the present work these two correlations were calculated from 
cross-correlations defined as: 
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( 11 ) 


(p';u')(x) = ^Jp'(t)u'(t + T)dt 

1 o 

1 T 

and, (p';u'u')(x) = — [ p'(t) [u'(t + x)] 2 dt 

T j 

0 

where x is the correlation time or delay time and T is the duration of the measured data. Note that by 
definition: 


p'u' = (p';u') T=0 and, p'u'u' = (p';u'u') T=0 


( 12 ) 


Fundamentals of Flow Measurement Using Rayleigh Scattering Principle 

A simplistic description of the measurement process, using laser induced Rayleigh scattering is 
schematically shown in figure 1 . When a laser beam is allowed to pass through a gas, the molecules present 
in the gas cause light scattering. The Rayleigh scattering process describes most (~99 percent) of the 
molecular scattered light. In the present experiment scattered light is collected and spectrally resolved to 
measure velocity. The line width of the incident laser beam is very narrow while that of the Rayleigh 
scattered light is wide. Even if the gas medium is stationary, the random thermal motion of the gas 
molecules creates a wide range of Doppler shift — resulting in a spectral broadening of the collected light. 
The Full Width at Half Maximum (FWHM) is a function of the distribution of molecular velocities and, 
therefore, is a measure of gas temperature. In the case of a moving gas media, the bulk velocity is 
superimposed on the random thermal motion of the individual molecules; therefore separation between 
the peaks of the incident laser line and the Rayleigh spectrum provides a measure of the bulk velocity. 

The total light under the Rayleigh spectrum is proportional to the molecular number density and provides 
a measure of gas density. Thus a single Rayleigh spectrum carries information of one component of gas 
bulk velocity, temperature and density. 

This basic principle has been used in the past to measure time-averaged quantities. Since density 
variation modulates the total scattered light, unsteady density fluctuations are easier to measure. 

Extension of the Rayleigh scattering technique to measure unsteady velocity fluctuation has remained a 
challenge. Recently, Seasholtz, Panda, and Elam (2002) have simultaneously measured the time variation 
of density and velocity fluctuations in free jets. The present setup is for a point-measuring system; a 
continuous wave laser was used, and scattered light from a small region was spectrally analyzed using a 
Fabry-Perot interferometer. 

To illustrate the unsteady velocity measurement process, first the nature of a fringe formed by the 
interferometer is shown in figure 2. The interferometer basically images the fiber-end delivering the 
scattered light. Since the fiber diameter is small, the field of view in the image covers a fraction of the 
free spectral range. The spectral analysis is a two-step process. First, a small portion of the light directly 
out of the laser beam is imaged through the interferometer (the Rayleigh scattered light is blocked; instead 
a part of the unscattered incident laser beam is analyzed). The narrow line width of the laser makes a 
sharp, narrow ring in this image (reference image, fig. 2(a)). In the second step the Rayleigh scattered 
light from a moving gas medium is passed through the fiber. The image formed by the Fabry-Perot 
interferometer (Rayleigh image, fig. 2(b)) is different: a diffused ring with a different ring diameter 
results. The radial shift in the peak intensity locations between the Rayleigh and the reference images is 
related to the Doppler shift associated with the bulk motion of the air stream (the diffused nature of the 
Rayleigh image is due to the thermal broadening). Since the laser frequency, and therefore, the reference 
fringe are fixed, an instantaneous change of air velocity produces a radial shift of the ring seen in the 
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Rayleigh image. To measure the time variation of velocity one needs to monitor the ring diameter in the 
Rayleigh image. To this end Seasholtz et al. (2002) used an image dissector that split the Rayleigh image 
into two concentric parts (fig. 2(c)), and measured the ratio of light intensities from the inner and the outer 
parts using two photomultiplier tubes (PMT). The ratio of photoelectron counts from the two PMT carried 
information about the ring diameter variation in the Rayleigh image, which in turn reflected the change of 
Doppler shift frequency from velocity change at the probe volume. For calibration, the intensity variations 
from the two parts were measured from known velocity flows. It was found that the intensity variations in 
either part of the image can be modeled by second order polynomials. If N 2 and N s denote photo-electron 
count rates from the inner and outer PMTs then: 


N 2 -A; +uB; +U 2 C-j , 

N3 =A 0 +uB 0 +U“ C 0 

Where, Aj, B,, C,, A 0 , B 0 and C 0 are calibration constants. The velocity component u is measured from a 
ratio of the two counts, R=N3/N 2 : 


~(RBj -B 0 )+^(RBj -B 0 ) 2 -4(RAj -A 0 )(RCj -C 0 ) 
2(RQ-C 0 ) 


(14) 


Note that the ratio of two counts cancels out changes in overall scattering intensity associated with flow 
density variation. Also the physically meaningful, positive root of the quadratic equation is considered. 
The component of velocity measured using a given optical arrangement depends on the angular position 
of the collection lens with respect to the incident laser beam. The present setup, shown in figure 3, 
measured the axial velocity component. Note that the setup measures negative -u velocity. The 
calibration process accounted for the sign reversal. The Rayleigh scattered light was collected and 
transmitted to an adjoining room for spectral analysis. Figure 4 shows a schematic of the spectroscopic 
arrangement. 

An important caveat in the above analysis is an implicit assumption that the effect of temperature 
fluctuations is also accounted for via the calibration process. Temperature broadening results in changes 
of light intensity in both the inner and outer parts of the image. By virtue of measuring a ratio of 
intensities between the two parts of the image, the process was made somewhat insensitive to temperature 
fluctuations. A numerical uncertainty analysis (Seasholtz et al. 2001) demonstrates that the effect of 
temperature fluctuation is small compared to the bigger change associated with velocity fluctuations. In 
addition, the present experiment was conducted in the potential core of the same free jets used to obtain 
calibration constants. Hence, the influence of isentropic temperature fluctuations was automatically 
included in the calibration constant. However, non-isentropic temperature fluctuations have affected the 
velocity fluctuations data in a small yet unknown way. 

To measure density fluctuations about 10 percent of the Rayleigh scattered light was split off using a 
beam splitter (BS1 in fig. 4) and the intensity fluctuations were monitored using a separate photo- 
multiplier tube. For a fixed optical setup and a fixed gas composition, the intensity variation of Rayleigh 
scattered light is directly proportional to air density variation. Once again, a calibration is necessary to 
relate photo-electron count to density change. The calibration was performed in the same core regions of 
the jet flows used for velocity calibration. If the local air density is p and photo-electron count rate from 
PMT1 is Ni then: 


Ni — C p i+ pCp2 • 


(15) 
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Where, C p[ and C p2 are calibration constants. In summary, a total of three PMTs were: two for measuring 
velocity and one for density. Additionally two sets of calibration constants were needed for velocity and 
density measurements. 


II. Experimental Set-Up 

Experiments were performed at NASA Glenn Research Center using three different nozzles (one 
convergent and two convergent- divergent) operated at Mach numbers, M = 0.95, 1.4 and 1.8. The 
convergent-divergent nozzles were designed by the method of characteristics. All nozzles were 25.4 mm 
in exit diameter. The jet facility used a continuous supply of unheated compressed air. The Rayleigh 
scattering system is somewhat elaborate and the following provides a brief description of some of the 
central features. An in-depth discussion of the Rayleigh set-up can be found in Seasholtz, Panda, and 
Elam (2002). The optical system was built in two parts. The first one is around the jet facility for 
transmitting laser light and collecting the scattered light (fig. 3). The scattered light was then passed 
through a 0.55 mm diameter optical fiber to a quiet room where the second part, consisting of a 
spectroscopic system and photon counting electronics, were placed (fig. 4). The splitting of the setup is 
necessary to minimize the effect of vibration on the optical components. Additional care had to be taken 
to reduce dust particles in air streams, and to stabilize the interferometer from temperature and vibration 
induced drifts. To reduce dust particles the dry air, supplied from a central high-pressure facility, was 
passed through additional micron size filters which made the primary jet air very clean. To clean the 
entrained ambient air, an additional air blower and filter system was installed that provided a large, 

200 mm diameter, low-speed (~20 m/s) co-flow around the 25.4 mm primary jet. 

The optical system was built over an X-Y traversing unit that carried laser head, transmission and 
collection optics. Surveys were made by moving the probe volume from point to point in the plume. The 
laser head of the solid-state, frequency-doubled, Nd:V0 4 laser was placed at the bottom part of the set-up. 
About 5 watts of single mode, 532 nm wavelength laser light was transmitted through a hollow side beam 
that contained a half-wave plate, focusing lens, mirrors and baffles. Since Rayleigh scattered light is 
polarization dependent, the half-wave plate was rotated to align the peak scattering plane with receiving 
optics. The background scattered laser light was significantly attenuated by suitable use of baffles and 
beam-dump. It was found that the noise from the jet created a tonal excitement of the laser line at around 
430 Hz. To reduce this excitation, an anechoic box was built around the laser head. This box significantly 
reduced the laser unsteadiness, but a trace remained and was manifested in the experimental results. 

To measure axial component of velocity the laser beam was passed at 45° to the jet axis, and scattered 
light from 90° to the incident direction was collected by lenses placed in the same plane containing the 
laser beam and the jet axis. The collection optics was made of a collimator (300 mm focal length, 82 mm 
diameter achromat) and a focusing lens (160 mm focal length achromat) that imaged the probe volume on 
the face of a 0.55 mm diameter multimode fiber. The combination of the fiber diameter and the 
magnification ratio of the collection optics defined the probe volume length to 1.03 mm. Not shown in 
figure 3 is an additional part of the set-up where a small part of light from the transmitting laser beam was 
split off for the purpose of monitoring laser frequency as well as to maintain the alignment in Fabry-Perot 
interferometer. Whenever necessary, a pneumatically actuated mirror was placed in the transmission laser 
path to divert light towards a diffuser. The mirror also blocked the incident laser beam from going 
through the jet, and the diffuser scattered the incident light; a small part of which was then collected by 
the optical fiber. 

In the second part of the setup (fig. 4) light arriving via the optical fiber was collimated by 1 00 mm 
focal length lens LI, and about 10 percent was split by BS1 and measured by PMT1. Output from PMT1 
was proportional to the air density fluctuations. The rest of the collimated beam was passed through a 
70 mm aperture planer mirror Fabry-Perot interferometer for spectral analysis. Single wavelength light 
from an extended source is imaged as concentric rings (fringes) at the output of the interferometer. 
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However, restriction of the field of view, imposed by the fiber diameter, created only one fringe as shown 
earlier in figure 2. The fringe-forming lens, which ultimately images the fiber face on the image dissector, 
was made of two camera lenses with suitable separation for an effective focal length of 2909 mm. The 
large magnification ratio of the setup created a 16 mm diameter image of the 0.55 mm diameter fiber on 
an image dissector. The image dissector was made of two concentric and slightly tilted mirrors. The inner 
one had a diameter of 1 0 mm and directed the inner part of the fringe to PMT2, while the 25 mm diameter 
outer one directed the outer part of the fringe to PMT3. The ratio of light intensities from PMT2 and 
PMT3 provided a measure of either axial or radial velocity as described earlier. 

The success of the velocity measurement system was critically dependent on stable operation of the 
interferometer. Slight thermal drift or change in the incident laser frequency displaced the reference fringe 
(fig. 2(a)), which was manifested as an artificial bias in velocity measurement. This made the auto- 
alignment setup (fig. 4) necessary. The alignment system was a feedback control that first measured the 
reference fringe diameter and compared it with a prescribed targeted diameter. This was accomplished by 
splitting parts of the transmitted light using a 3-prism assembly, and imaging them on a CCD camera. 
Subsequently, the difference between the targeted and the measured fringe diameters was translated into a 
change of high voltage supply to the piezo-electric actuators that adjusted interferometer’s plate 
separation. Before every Rayleigh measurement, reference light was collected and the auto-alignment 
system was engaged. When the desired fringe diameter was obtained within a tolerance, the reference 
light collection system was disengaged, and velocity and density measurement via analysis of Rayleigh 
scattered light began. 

Photoelectron counting electronics were used with all 3 PMT signals. The photoelectron pulses from 
each PMT arriving within a preset gated time interval was counted using a PC counter-timer board. 
Counting was performed on contiguous series of gated time intervals without any dead time between the 
gates. Long time records of up to 5 million data points from each PMT were acquired with a typical 
sampling rate of 90,000/sec. The data points were converted to physical parameters: instantaneous density 
and velocity through the usage of proper calibration constants. The Welch method of modified 
Periodograms (Welch, 1967) was used to calculate individual power spectra and cross-spectra. Each long 
record was divided into small, 50 percent overlapped, segments; modified periodograms of each segment 
provided local estimates; average of all local estimates provided the final power spectra. Time domain 
correlation calculations were also performed via Fourier transform, where individual segments of velocity 
and density time histories were Fourier transformed, multiplied and inverse transformed. For example, the 
density-velocity cross-correlation was calculated as: 

(p';u') = 3- 1 (3(p(t))*3*(u(t))) (16) 

In this equation 3 represents Fourier transform, complex conjugation is represented by the superscript *, 
and inverse Fourier transformation by 3 _1 . Average of local estimates produced the final correlation. The 
length of the data segment was varied to assure independency of the calculated values on string length. 


III. Results 

Validation of the Technique 

Figure 5 presents a comparison of the velocity fluctuations spectra measured using the Rayleigh 
technique and a hot-wire probe. This figure shows that the shapes of velocity spectra are similar while the 
absolute energy level in the spectrum measured by the Rayleigh technique is two and a half times of that 
measured by hot-wire. The Rayleigh spectra float on constant white noise floors, which are expected 
consequences of electronic shot noise. Electronic shot noise creates the unavoidable, fundamental 
uncertainty in all optical intensity measurements. The photoelectron count in every time bin is 


NASA/TM— 2005-213564 


7 



contaminated by random fluctuations that follow Poisson distribution. Spectra calculated from the time 
series of photoelectron counts obtained from still air, where there is no density or velocity fluctuations, 
show a white noise floor. The broadband, frequency independent, flat spectrum is a property of electronic 
shot noise. The propagation of electronic shot noise in the velocity spectrum is somewhat complicated by 
the ratio of photoelectron counts, R = N 3 /N 2 , needed to determine instantaneous velocity. Nevertheless, 
the above spectral property was used to extract standard deviation of the velocity fluctuations from the 
spectra. The measured fluctuations spectra S u -2 is due to a sum of the contributions from the shot noise 

(with mean square value O“ hot ) and the desired velocity fluctuations (mean-square uj-ms)- The total 
energy under the power spectrum can be expressed as (Saleh and Teich 1991): 

Js u . 2 df=u& s+0s \ ot (17) 

The constant shot noise floor was estimated as the spectral density value at the highest resolved frequency 
(Nyquest limit) of 45KHz: S s | 10t = (S u -2 ) f=4500Q where it was known that the energy from turbulent 

fluctuations was small. Figure 5 shows the estimated noise floor for each Rayleigh spectra. 

SinCea shot = j ( S U ' 2 )f =45000 df ’ ( 18 > 

the standard deviation of velocity fluctuations is: 

U nns = -\Zj[S u '2 ~(S U '2 )f =45000 ]df • 

The rms values measured by this shot noise subtraction process are found to be reasonable. For example, 
in figure 7(a) (u / rm s/Uj) h ot -wire 0.166 and (u rm s/Uj)R a yleigh corrected 0.154 and in figure 7(b) (u^/UjVt 

-wire 

0.14 and (Um/UjiRayieigh conected = 0.132. The subtraction method was used uniformly to estimate p nns and 
u'nns from respective spectra. 

Figures 6(a), 6(b), 7(a), and 7(b) show typical velocity and density fluctuations spectra obtained from 
a region of very high turbulent fluctuations in Mach 1.4 and 0.95 plumes. Flot-wire is unusable in these 
jets and spectra measured by Rayleigh technique are shown. It is to be pointed out that the presented 
velocity and density spectra were narrow-band filter around 430 Flz. The unfiltered velocity data show a 
large spike at this frequency. The spike was found to be independent of the jet speed and was attributed to 
a spurious oscillation in the incident laser frequency. As mentioned earlier, the laser head was placed 
directly underneath the jet and the high noise level created by the plume is suspected to cause small 
vibration in the laser cavity leading to the spike (Panda et al. 2003). Nonetheless, the density and velocity 
fluctuations spectra are similar in shape. For data in figures 6 and 7 the spectra show significant energy in 
the lowest frequency. The high frequency roll-off is masked by the shot noise floor. 

The cross-correlation calculations between density and velocity fluctuations, however; are mostly 
unaffected by shot noise. Since the shot noise from different PMT are uncorrelated, the cross-correlation 
process eliminates the influence of electronic shot noise and the absolute cross-correlation values are 
expected to have very small uncertainty. Parts (c) and (d) of figures 6 and 7 presents the non- 
dimensionalized cross-correlation values measured in two different Mach number jets. Note that the 
correlations are expected to be symmetric in time; only the positive time delays are shown here. There are 
multiple interesting observations that can be made from these figures. First, <p ;u> correlations are 
always fairly high and mostly positive. Second, <p / ;uu> correlations can be high or low depending on the 
Mach number and probe locations. In the present unheated jets positive values of the <p ;u > correlation 
are expected. The unheated jets cool to lower temperatures and higher densities than the ambient upon 
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expansion through the nozzle. Since the high velocity is associated with the high density (jet air) and the 
low velocity with the low density (entrained ambient air) <p ;u 7 > correlations are positive. It is expected 
that if the same experiments were conducted in a heated jet, with core density lower than the ambient, 
negative correlations would have resulted. The <p ;u u > correlations, on the other hand; do not follow this 
argument. Following is a discussion of the uncertainty in the cross-correlation data. 

The correlation values are normalized by individual standard deviations which are the primary 
sources of uncertainty in the presented form. For example, correlations between velocity and density 

(p';u') , , 

fluctuations are expressed as — y y — , which requires estimates of u ms and p ms . As mentioned above 

Prms u nns 

<p;u> correlations are relatively error free, however; the shot noise subtraction process described earlier 
has to be used to determine u mi; , and p rnK . The subtraction process creates ±10 percent error in the 
standard deviation estimates. Flence, uncertainty levels of ±20 percent are expected on the quoted 
normalized correlation values. This large uncertainty supersedes all other sources, such as the 
contamination from temperature fluctuations, numerical convergence error etc. 

The differences between Favre averaged and Reynolds averaged values are straightforward to 
calculate from the measured correlation data. For example, for the conditions shown in figure 7 the 
normalized differences between the Favre-averaged and time-averaged axial velocity is calculated from 
equations (5) and (12), using additional measured values of p = 1.44 Kg/m 3 , u = 391m/s and 

u' 2 = 1028 (m/s) 2 : 


u-u _ ( p ,u )t= 0 
u pu 


0.004. 


Note that the large uncertainty in the normalized <p / ;u / >/(p / rms u / rm s) does not affect the above estimate. 
Since the absolute values of the correlations are measured relatively accurately, the primary source of 
error is the uncertainty in p and u measurement. It is estimated that (u - u)/u is determinable within a 

maximum of ±10 percent uncertainty. The difference in the kinetic energy associated with the axial 
velocity fluctuations is calculated using equations (10) and (12): 


uV'-u'u' 1 

(p" u ' u ') T =o 

f»Ao] 

2 n 

u'u' u'2 

p 

P 

V J 

/ 


= -0.038 


The uncertainty in the above calculation is relatively higher due to the uncertainty in the u' 2 estimates. 
The estimated level is ±25 percent of the quoted number. 

A large number of data points were measured from many axial and radial locations in different Mach 
number plumes. Data analyzed from some selected conditions are shown in table 1. A closer look into this 
table confirms previous observations on <p ;u > and <p / ;u / u / > correlations. Additionally it confirms that 
the Favre averaged velocity is very slightly (fraction of a percent to one percent) higher than the Reynolds 
averaged counterpart. Favre-averaging also produces lower level (fraction of a percent to four percent) of 
kinetic energy compared to the Reynolds averaging. Notably <p ;u u > correlations have stronger impact 
than <p ;u y > correlations on the differences in the kinetic energy; a larger value of the former reflects as a 
larger differences between the two averaging process. 
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Table 1 also provides a comparison of the measured air density fluctuations p rms / p with that 
estimated from velocity fluctuations u^u using Morkovin’s hypothesis (eq. (3)). Morkovin’s hypothesis 
is always found to under-predict fluctuations at all Mach number conditions. The differences are 
particularly large in regions with vigorous turbulent fluctuations such as the end of the potential core and 
the middle of the lip shear layer. 


IV. Summary and Conclusion 

A molecular Rayleigh scattering based optical measurement technique was used to simultaneously 
measure air density and the axial component of velocity at a sampling rate of 90,000/sec. The density 
variation was measured by monitoring the variation of the total scattered light intensity. The velocity 
variation was measured by spectrally resolving the scattered light using a Fabry-Perot interferometer and 
monitoring the time dependent Doppler shift from the incident laser frequency. A small probe volume 
length of 1.06 mm and diameter 0.2 mm 2 allowed point measurement from different regions of the flow- 
field. One convergent and two convergent-divergent nozzles were used to produce unheated plumes in 
Mach numbers 0.8, 0.95, 1.4 and 1.8. The standard deviations of the density p^ and axial velocity u rms 
fluctuations and the density-velocity <p ;u 7 > and density-velocity 2 <p / ;u u / > cross-correlations measured 
from different regions of the plumes provided the following results. 

(a) Normalized p'u' correlations are significantly high (varies between 0.2 and 0.7) and almost always 
positive. The positive correlation is expected in the unheated plume where the jet core has higher density 
than the ambient. 

(b) Normalized p'u'u' correlations are also high (mostly varies between -0. 1 and -0.5) and always 
negative. 

(c) The Favre averaged axial velocity is very slightly (fraction of a percent to one percent) higher than 
the Reynolds averaged counterpart. 

(d) The axial component of turbulent kinetic energy calculated via Favre-averaging, 0.5 u "u " , is 

lower by a fraction of a percent to four percent from the same calculated via Reynolds averaging 0.5u'u' . 

(e) Estimation of turbulent density fluctuations, from the axial velocity fluctuations using Morkovin’s 
hypothesis, is found to provide lower values compared to the measured data at all Mach number 
conditions. 

The present experiments were conducted in unheated jets where the stagnation temperature of the 
plume was the same as that of the ambient. It will be of interest to conduct similar experiments in heated 
jets where the effect of total temperature fluctuations on the above parameters can be established. 
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Figure 1. — Schematic of Rayleigh spectrum from air flow. 


Images of the fiber face 



0^ 


( 

o 

i 



w 
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by PMT2 
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Figure 2. — Fringes formed after passing through Fabry-Perot Interferometer by (a) incident 
laser light (b) Rayleigh scattered light; (c) splitting of Rayleigh 
image to measure velocity fluctuations. 
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Fiber to loft 



(a) 



(b) 


Figure 3. — (a) Top view of optical arrangement to measure 
axial u component of velocity, (b) Side view. 


Auto-alignment 



Figure 4. — Schematic of optical setup to analyze collected light. LI to L5 are lenses; 
BS1 is beam-splitter and FPI is the Fabry-Perot interferometer. 
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Figure 5. — A comparison of power spectral density of u fluctuations 
measured by the Rayleigh technique and a hot-wire 
probe at indicated locations of Mach 0.6 plume. 
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Figure 6. — Power spectral density of (a) air density fluctuations and (b) axial velocity 
fluctuations. Correlation between (c) density-velocity fluctuations, 
and (d) density-velocity 2 fluctuations measured 
at x/D =10, r/D=0 in a Mach 1 .4 jet. 
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Figure 7. — Same as figure 6 except M=0.95 jet, probe 
volume location x/D=10 and centerline. 
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